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Abstract — The  performance  of  the  conventional  bearings-only 
tracking  (BOT)  from  a  single  passive  sensor  largely  depends  on 
the  sensor  platform  maneuvers.  This  paper  presents  a  new  BOT 
approach  based  on  fusion  from  two  heterogenous  bearings-only 
sensors  residing  on  the  same  moving  or  stationary  platform.  The 
two  sensors  are  an  ESM/EO  with  negligible  propagation  delay 
and  an  acoustic  sensor  with  significant  propagation  delay.  Since 
target  range  information  is  contained  in  the  acoustic  propagation 
delay,  the  problem  is  therefore  observable  even  when  the  platform 
is  stationary.  An  out-of-sequence  measurement  fusion  from  the 
acoustic  and  ESM/EO  sensors  (OOSM-AE)  is  developed  to 
estimate  the  target  trajectory.  It  consists  of  an  unscented  Kalman 
filter  (UKF)  to  handle  in-sequence  ESM/EO  measurements  and 
an  OOSM  unscented  Gauss-Helmert  filter  (OOSM-UGHF)  to 
handle  out-of-sequence  acoustic  measurements.  Simulation  tests 
are  conducted  to  demonstrate  the  performance  of  this  new  BOT 
approach. 

Keywords — Propagation  delay,  bearings-only  tracking,  target 
motion  analysis,  unscented  Gauss-Helmert  filter,  out-of-sequence 
measurement. 

I.  Introduction 

The  commonly  used  passive  sensors,  like  acoustic  sensors, 
electronic  support  measures  (ESM)  sensors  and  electro-optical 
(EO)  sensors,  detect  target  bearings  only.  This  makes  the 
target  trajectory  estimation  from  range-absent  measurements 
a  challenging  problem. 

Several  approaches  have  been  developed  in  the  last  four 
decades.  The  most  popular  one  is  to  deploy  a  passive  sensor  on 
a  maneuvering  platform,  and  the  target  trajectory  is  estimated 
using  bearings-only  tracking  (BOT)  or  bearings-only  target 
motion  analysis  (BO-TMA)  [15]  [1].  This  approach  requires 
the  sensor  platform  to  maneuver,  so  the  target  trajectory  is 
observable  [16]  [11]  [6].  Since  these  maneuvers  can  interfere 
with  the  sensor  platform’s  own  mission  (for  example:  to  reach 
its  destination  as  early  as  possible),  BOT  from  a  nonmaneuver¬ 
ing  sensor  has  attracted  attention  recently.  Results  showed  that 
the  BOT  problem  is  indeed  observable  from  a  nonmaneuvering 
sensor  when  a  target  is  performing  particular  maneuvers  (two- 
leg  with  constant  speed,  or  constant  turn)  [12]  [7].  However, 
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there  is  still  a  gap  to  transition  these  results  to  real  applications, 
for  the  target  can  maneuver  in  a  manner  unbeknownst  to  the 
observer,  and  the  nearly  constant  velocity  (CV)  motion  is  the 
common  motion  used  in  operation. 

The  BOT  approach  has  been  extended  to  the  Doppler- 
bearing  tracking  (DBT)  approach  in  [17]  [10],  This  approach 
tracks  target  position  and  emitted  frequencies  from  bearings 
and  Doppler  shifted  frequencies  and  a  target  trajectory  can 
be  estimated  even  when  the  platform  is  not  maneuvering. 
The  difficulty  faced  in  DBT  is  to  identify  target  frequencies 
from  a  noisy  environment,  especially  when  the  target  emitted 
frequencies  are  varying. 

Another  approach  is  to  locate  targets  through  triangulation 
from  multiple  stationary  or  moving  passive  sensors  located  at 
different  positions.  This  approach  needs  to  remove  triangula¬ 
tion  ’’ghosts”  in  multi-target  scenarios,  and  can  be  solved  as  an 
S-D  assignment  problem,  where  S  is  the  number  of  sensors. 
A  Lagrangian  relaxation  was  suggested  to  solve  the  problem 
when  S'  >  3  [18]  [8].  By  making  use  of  Doppler  frequencies, 
the  number  of  sensors  can  be  reduced  to  2  (S  =  2)  [20]. 

In  this  paper,  we  propose  a  new  bearings-only  approach 
to  fuse  measurements  from  two  heterogenous  passive  sensors 
deployed  on  the  same  platform  which  can  be  either  moving 
or  stationary.  The  two  sensors  are  a  passive  ESM/EO  sensor, 
designated  as  sj  and  a  passive  acoustic  sensor,  designated  as 
S2-  Both  sensors  detect  target  bearings  only.  The  ESM/EO 
sensor’s  detections  have  no  propagation  delay,  whereas  the 
acoustic  sensor  receives  the  target  signals  after  significant 
propagation  delays.  The  time  difference  between  the  reception 
times  of  the  two  sensors  is  the  acoustic  propagation  delay,  and 
the  target  range  can  then  be  inferred  from  the  estimates  of  these 
delays  assuming  the  propagation  speed  is  known.  Complete 
observability  in  this  BOT  problem  is  therefore  obtained,  as 
range  is  implied  in  the  sensors’  reception  times. 

However,  to  obtain  target  range  using  the  principle  men¬ 
tioned  above  is  not  straightforward.  To  compute  the  acoustic 
propagation  delay,  a  pair  of  passive  signals  from  si  and  S2 
having  the  same  emission  time  needs  to  be  identified.  A  BOT 
target  usually  emits  continuous  signals  which  are  received 
by  the  sensors  and  discretized  by  sampling.  They  are  not 
instantaneous  signals,  like  ’’ping”  or  ’’pulse”  which  can  be 
associated  easily.  There  is  no  feature  to  identify  an  acoustic 
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bearing  measurement  and  an  ESM/EO  measurement  emitted  at 
the  same  time.  Furthermore,  the  ESM/EO  and  acoustic  sensor 
may  have  different  sampling  times  (they  are  asynchronous), 
and  the  sensor  platform  may  be  dynamic.  These  make  the 
problem  even  more  complicated.  A  comprehensive  algorithm 
is  therefore  needed  to  take  all  these  issues  into  consideration. 


Recently,  we  have  formulated  implicit  constraint  dynamic 
estimation  problem  to  a  Gauss-Helmert  model  (GHM),  and 
proposed  an  unscented  Gauss-Helmert  filter  (UGHF)  [22] 
[21]  to  solve  this  problem.  The  existing  UGHF  works  with 
in-sequence  measurements.  Further  development  on  OOSM- 
UGHF  is  presented  in  this  paper. 


Acoustic  is 2) 


ESM/EO  (i,) 


t  1  t  1 
lk- 4  lk- 3 


I  ► 

t‘l  ]  reception  time 


Fig.  1.  Out-of-sequence  measurements  in  ESM/EO  and  acoustic  sensors. 


Fig.  1  illustrates  the  ESM/EO  and  acoustic  signal  emission 
and  reception  time  sequences,  where  k  is  the  reception  time 
index,  which  orders  the  combined  acoustic  and  EO/ESM 
discretized  signals  by  arrival  (sensor)  time.  This  is  also  the 
measurement  index,  while  i  and  j  are  the  target  signal  emission 
time  indexes  from  si  and  s2,  respectively.  It  can  be  seen 
that  out-of-sequence  measurements  (OOSM)  occur  due  to  the 
acoustic  propagation  delay. 

The  OOSM  problem  is  also  referred  to  as  “negative¬ 
time  measurement  update”  problem,  namely,  the  state  time, 
i®2,  corresponding  to  the  latest  measurement  at  tsk 2  is  earlier 
than  the  latest  state  updating  time,  t®1,  namely  t®2  <  f®1 
The  prediction  step  in  the  in-sequence  estimation  becomes  a 
retrodiction  in  the  OOSM.  The  OOSM  problem  has  been  well 
studied  [5].  The  simplest  approach  performs  an  approximate 
retrodiction  by  neglecting  the  process  noise  [5].  This  approach 
is  referred  to  as  Algorithm  C  [5].  Algorithms  B1  and  A1 
were  proposed  to  solve  the  one-step-lag  OOSM  by  considering 
process  noise  [9]  [3],  and  they  give  an  approximate  and  the 
exact  solutions,  respectively.  They  were  further  developed  to 
the  algorithms  B/l  and  A/1  for  solving  the  /-step-lag  OOSM 
(/  >  1)  in  a  single  step  [4]. 

The  existing  OOSM  algorithms  mentioned  above  assume 
that  retrodiction  time  is  known.  However,  the  retrodiction  time 
is  the  acoustic  signal  emission  time  in  our  problem.  This  is 
unknown  to  the  observer  and  depends  on  the  state  of  the  target 
which  is  given  by  the  following  propagation  delay  constraint: 

tr=tsi?-s>  a) 

where 

Si  =  %  (2) 

is  the  propagation  delay,  cp  is  the  signal  propagation  speed 
in  the  medium,  and  rj,  which  depends  on  the  state  (at  the 
emission  time),  is  the  distance  from  the  target  at  time  i®2  to 
the  sensor  at  time  tsk2 .  This  leads  to  an  implicit  constraint  in 
state  transition  model. 


The  aim  of  this  paper  is  to  develop  a  comprehensive  al¬ 
gorithm  to  estimate  states  with  fusion  of  in-sequence  bearings 
from  the  ESM/EO  sensor  (si)  and  out-of-sequence  bearings 
from  the  acoustic  sensor  (s2).  State  estimation  for  the  bearings 
from  si  will  be  performed  by  a  unscented  Kalman  filter  (UKF), 
and  for  s2,  a  new  OOSM-UGHF  will  be  developed.  They  are 
described  in  the  following  sections. 

Unscented  transform  (UT)  is  selected  to  approximate  non¬ 
linear  transformations  in  both  in-sequence  and  out-of-sequence 
estimations  in  this  paper.  This  is  because  it  provides  better 
accuracy  than  the  first-order  Taylor  linearization  (used  in  the 
extended  Kalman  filter)  by  accounting  for  the  asymmetry 
of  the  nonlinear  transformation.  Compared  to  the  particle 
filtering  approach,  a  UT  is  orders  of  magnitude  less  demanding 
computationally.  The  UT  has  advantages  in  both  estimation 
accuracy  and  efficiency,  which  are  important  factors  in  real 
applications. 


II.  State  Estimation  with  the  ESM/EO  Sensor 


The  state  estimation  for  the  ESM/EO  bearings  is  straight¬ 
forward  as  the  measurement  is  in-sequence  and  no  propagation 
delay  needs  to  be  taken  into  consideration.  The  problem  is 
formulated  based  on  the  nearly  CV  state  model  (or  WNA  — 
white  noise  acceleration).  The  target  state  with  size  4  is  defined 
as 

=  [  *(#)  !/(#)  *(#)  »(#)]'  O) 


where  tsk  is  the  signal  reception  (or  sensor)  time  by  the 
ESM/EO  sensor  si  at  time  cycle  k.  Since  the  propagation  delay 
is  negligible  for  si,  the  target  signal  emission  time  t'-1  is  equal 
to  tsk  .  The  state  transition  model  is1 

x4^1) = m1  ^-1)^-1) +v4(«i-1)  (4) 


where  the  transition  matrix  is 
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with  v4  the  zero-mean  process  noise  (WNA)  for  the  interval 
(Tfc1 .  f^.1— 1].  The  resulting  discretized  white  noise  acceleration 
(DWNA)  model  [2]  has  covariance 
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1  f  fere  it  is  assumed  for  simplicity  that  the  measurements  arriving  at  tk.  | 
and  are  both  from  sensor  s\. 
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where  q  is  power  spectral  density  (PSD)  of  the  motion  process 
noise  (same  for  x  and  y,  and  assumed  independent  between 
the  coordinates).  The  measurement  model  is  given  by 


z(tk)  =  Htk)=t  an  1 


Uo-m). 


+w(tsk ) 


(8) 


where  xs(tsk)  and  ys(t] ())  are  the  sensor  positions  at  time 
tsk  in  the  x  and  y  coordinates  respectively,  w (tsk)  is  zero- 
mean  white  Gaussian  measurement  noise  with  variance  R(tk' ), 
assumed  independent  of  the  process  noise. 


The  unscented  Kalman  filter  (UKF)  is  used  to  estimate  the 
state  [13]. 


III.  State  Estimation  with  the  Acoustic  Sensor 

An  out-of-sequence-measurement  filter  is  required  for  the 
bearings  from  the  acoustic  sensor  s2-  It  can  be  seen  in  Fig.  1 
that  an  acoustic  measurement  received  at  time  tsk  corresponds 
to  the  target  state  at  emission  time  £®2,  which  is  earlier  than 
the  latest  state  assumed  to  have  been  updated  by  the  ESM/EO 
sensor  at  time  tsk_x  =  t®1.  The  problem  is  then  to  update  the 
state  estimate  x4(ffc!_i|ffc!_i)  with  the  acoustic  measurement 
z(tk2).  The  main  challenge  of  this  problem  compared  to  the 
existing  OOSM  approaches  is  that  the  time  te2  is  unknown, 
and  it  needs  to  be  estimated  together  with  the  kinematic  state. 

The  OOSM  approach  to  address  the  above  mentioned 
problem  consists  of  the  following  steps: 

•  Retrodict  the  state  from  time  t®1  =  tsk1_1  to  the 

(unknown)  emission  time  f®2  (to  be  estimated)  cor¬ 
responding  to  the  sensor  time  tsk2 .  The  state  before 
retrodiction  is  and  the  state  after  retro- 

diction  is  x5(f®2 The  latter,  defined  in  (9), 
includes  the  acoustic  emission  time. 

•  Update  the  state  estimate  x4(i£_i|tfc!_i)  to 
x^ffcLil  tsk)  by  the  acoustic  OOSM  z(tsk). 

The  unscented  transform  is  used  in  the  above  two  steps  instead 
of  the  first-order  Taylor  linearization  used  in  the  existing 
OOSM  algorithms  [4]  [5], 

A.  State  Retrodiction 

The  retrodiction  has  to  be  done  to  the  emission  time  te-2 
that  is  unknown  to  the  observer,  but  can  be  derived  from  the 
propagation  delay  constraint  described  in  (1).  To  estimate  the 
retrodicted  target  kinematic  information  and  the  emission  time 
f®2  simultaneously,  the  following  augmented  state  is  defined: 

x5(f  ?)=[x{t?)  y{tf)  x(t  f)  y(tf)  tf]'  (9) 

Obviously,  the  positions  x{te2),  y{tj2)  and  the  time  P2  de¬ 
pend  on  each  other,  and  this  leads  to  the  retrodicted  state 
x5 {tj2  \t3k_1)  and  the  latest  state  estimate  x4(£^1_1  I^Li )  to 
have  an  implicit  relationship.  The  Gauss-Helmert  transition 
model  [22]  [21],  which  handles  such  implicit  relationships, 
is  then  used  for  retrodiction.  This  is  described  by 

g  [xB(tf  J.X4^!)]  +vB(i£_1,i?)  =  05  (10) 

where  g[  ]  is  the  Gauss-Helmert  implicit  state  transition  func¬ 
tion,  which  combines  the  target  motion  constraints  and  the 


delay  constraint  between  x5(P2)  of  dimension  5  and  x4(if1_1 ) 

of  dimension  4.  Assuming  the  target  motion  follows 

a  WNA 

motion,  g[-j 

is  given  by 

g['] 

=  [  fli(')  52(0  gsi')  54(0  Ss(')  ]' 

(11) 
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and 

Tj,k-i  —  tj2  tk-i  <  0 

(17) 

r3 

=  yWf )  -  zs0fc2)]2  +  WT>  ~  ys(fk)}2 

(18) 

Note  that  Eq.  (16)  is  the  one  that  connects  the  emission  time 
and  location  to  the  corresponding  sensor  reception  time. 


The  process  noise  v5  is  the  zero-mean  Gaussian.  Based  on 
the  DWNA  model  [5],  its  covariance  is 


Q5(t?,tsk-i)  = 
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(19) 


where  q  is  as  in  (7),  and  qs  is  the  variance  of  the  process  noise 
in  the  delay. 

The  algorithm  used  for  retrodiction  is  the  UGHF  [22]  [21], 
which  obtains  the  retrodicted  state  iteratively  through  a  Gauss- 
Newton  algorithm.  Given  and  its  error  covari¬ 

ance  P4(£^1_1|£^1_i),  the  sigma  points  and  their  corresponding 
weights  are 

SigmaPts  [x4 (f x  |4L 1 ) ,  P4 (i£_  1  \t ) ,  «](20) 

where 
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where  m  =  —4, . . . ,  4,  is  the  sigma  point  index, 

L  / +  K)P{tl_i\tsk1_i)  \  indicates  the  |m|th  column  of 

4  +  /t)P 4(tfc1_1|f^1_1)  ],  and  k  is  a  scalar  that 
determines  the  spread  of  sigma  points.  Each  sigma  point  is 
retrodicted  from  the  previous  target  time  tsk1_l  to  an  unknown 
time  The  problem  is  then  to  solve 

g  |C_1),x4-m(^_1|^_1)]  =  05  rn~.  4, ....  4 

(26) 

Note  that  the  process  noise  is  not  taken  into  consideration  in 
the  OOSM  algorithm  C. 

A  Gauss-Newton  algorithm  is  applied  to  obtain  the  points 
*5,m(*TI*E-i)  iteratively-  The  iteration  procedure  (with  index 
n)  for  the  mth  sigma  point  is 

[*6,m(W-i)]B  =  l^-r)]”-1  +  A-1 

xg  [[&6”B(*5aiei)]n-1.*4,m(*i1-1IC1)]  (27) 

where  A  (without  arguments,  for  conciseness)  is  the  Jacobian 
matrix  given  by 
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The  initial  of  the  mth  sigma  point  [x 

5’m(fj2|ffc1_1)]°  for  the 

iteration  process  in  (27)  is  computed  approximately  by 


[xm(tf\tskU)]°  =  *m(C-i  l41-1)+im((fc1-1l41-1)x[A(fJe2)]0 

(32) 


(33) 

[xm(t?KU)]°  =  ^K-iK-d  (34) 

IlHW-r)]0  =  ym(tSk-iK-i)  (35) 

[{te2)m]0=tsk2_5j  (36) 

where 

[A(tf)]°  =  [(i?)TO]°-i?_i  (37) 

«  ^r/cP  (38) 


and  r^l_1  is  the  distance  between  the  target  estimate  and  the 
sensor  at  time  tsk1_1. 


The  retrodicted  state  x5(f^2|f^.1_1)  and  its  error  covariance 
P5(4j2|(fc-i)  are  then  computed  from  weighted  sums  of  the 
retrodicted  sigma  points. 
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771=  —  4 

puffier) 
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«  53  wmx5>m(f®2|f*i_1 

m=— 4 


(40) 

where 

*-5’m(tr K-i) = x5-™(ff i^_i)  _  x5(tf itjLj  (4D 
with  m  =  —4, . . . ,  4. 

B.  State  Update 

This  step  updates  x4(f^)_i|f^(_i)  to  x4(f^Lil^2)  by  the 
OOSM  z[tsk  )  —  it  fuses  the  latter  into  the  former.  Note  that 
the  sigma  points  of  x4(tj?Li|ffc)_i)  have  been  generated  in  (20). 

According  to  the  MMSE  estimator  [2],  x4(f^1_1|f^2)  and 
its  error  covariance  P4(tsk1_1\tsk)  are 

*4(*fc-i  1C)  =  +  p  xzP^zWk)  ~  ^k)] 

(42) 

P4(^-il42)  =  P  W-iK-i)  ~  P«P«P«  (43) 

The  expected  measurement  z(tsk)  is  based  on  the  retrodicted 
state  x5’m(fj2|ffc1_1)  as 
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The  variance  Pzz  of  the  innovation  and  the  covariance  Pxz 
are  computed  by 
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=  *m(^2)  -  ^2) 

(49) 

The  OOSM-UGHF  does  not  create  new  states,  it  only 
updates  the  states  generated  by  the  bearings  from  sj  before. 
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IV.  Simulation  Results 

Simulation  results  are  preseued  to  demonstrate  the  new 
algorithm’s  performance.  The  conventional  BOT  approach  is 
also  evaluated  using  the  same  simulation  data.  Two  sensor 
platform  scenarios  are  used  in  the  simulation  tests: 

•  Maneuvering  (M):  It  has  three  legs  linked  by  two 
90°  turns  with  turn  rate  3°/s  shown  in  Fig.  2.  The 
platform  speed  is  lOm/s  throughout  the  whole  path. 
It  moves  to  the  east  for  60s,  spends  30s  to  make  a 
90°  left  turn,  and  moves  to  the  north  for  60s.  It  then 
makes  a  90°  right  turn,  and  moves  towards  the  east 
for  180s.  The  total  duration  is  360s. 

•  Stationary  (S):  The  platform  stays  at  position  (0m, 
0m)  for  360s. 

An  ESM  sensor  and  an  acoustic  sensor  are  deployed  on  the 
platform  to  detect  target  bearings  regularly.  The  two  sensors 
are  not  synchronized.  Their  sampling  intervals  and  initial 
detection  times  are  different.  The  ESM  sensor  is  initiated 
at  time  0s  with  sampling  interval  Is,  whereas,  the  acoustic 
sensor  is  initiated  at  time  0.2s  with  sampling  interval  2s. 
The  measured  bearing  errors  of  the  ESM  and  the  acoustic 
sensors  are  zero-mean  white  Gaussian  with  standard  deviations 
a b  =  1°.  Both  sensors  have  no  bearing  detection  during  turns 
(total  missed  detection  duration  is  60s  for  the  maneuvering 
platform  scenarios). 


x  (km) 


Fig.  2.  Test  scenarios.  Initial  locations  of  the  targets  and  the  maneuvering 
sensor  platform  are  shown  as  ”o”.  The  stationary  platform  is  not  shown  in  the 
figure. 


both  stationary  and  moving  (maneuvering  or  nonma¬ 
neuvering)  platform. 

•  UKF-E:  A  UKF  to  estimate  state  based  on  the  ESM 
bearings  only.  The  acoustic  bearings  are  regarded  as 
“expired"  information  and  discarded.  This  algorithm 
is  considered  as  a  conventional  BOT  approach  which 
works  for  maneuvering  platform  only. 

The  initial  state  estimate  is2 


^(fg1)  =  [  f0sin&o  f0cosbo  0  0  f  (50) 


where  r0  is  set  to  10000m  which  is  more  than  4000m  away 
from  the  ground  truth  (ro  =  4000m~  6000m),  and  bo  =  b( fg1 ) 
is  the  ESM  measured  bearing  at  time  fg1  =  0s.  The  initial  state 
error  covariance  is  computed  by  [19] 


where 
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Pxx  =  {ro&b  cos  b0)2  +  (oy  sin  b0)2 

Pyy  =  (f0(Tb  sill  b0)2  +  (oy  cos  b0)2 

Pxy  =  Pyx  =  (of  -  r^al)  sin  b0  cos  b0 


(51) 


(52) 

(53) 

(54) 


where  ay  =  fo/3  is  the  initial  range  error  standard  deviation. 
The  process  noise  PSD  q  in  (7)  is  set  to  0.01m2/s3.  The 
acoustic  propagation  speed  cp  in  the  air  is  344m/s.  The  scalar 
k  in  (20)-(25)  is  set  to  1. 


Four  targets  moving  at  constant  speeds  of  5m/s,  lOm/s, 
50m/s  and  lOOm/s  respectively  are  shown  in  Fig.  2.  The  state 
estimation  starts  50s  later  after  targets  move  from  their  initial 
positions,  so  the  acoustic  signal  can  be  guaranteed  to  reach  the 
sensor  platform  when  estimation  starts.  This  means  that  the 
targets  are  at  their  initial  points  at  time  -50s,  and  the  sensor 
platform  is  at  its  initial  point  at  time  0s.  The  estimation  starts 
at  time  0s. 

The  algorithms  used  in  the  simulation  are: 

•  OOSM-AE:  The  acoustic-ESM  fusion  algorithm  pro¬ 
posed  in  this  paper.  The  OOSM-UGHF  is  applied  to 
the  bearings  from  the  acoustic  sensor,  and  the  UKF 
is  used  to  the  bearings  from  the  ESM.  It  works  for 


Fig.  3.  Maneuvering  platform:  The  estimated  position  RMSE  versus  time 
for  the  target  with  speed  5 m/s. 

The  simulation  results  were  obtained  from  100  Monte 
Carlo  runs.  The  estimated  position  root  mean  square  errors 
(RMSE)  versus  time  are  displayed  in  Figs.  3-8,  where  Figs.  3— 
6  are  for  the  maneuvering  platform,  and  Figs.  7-8  are  for 
the  stationary  platform.  The  overall  and  the  last  point  position 
RMSEs  for  all  the  scenarios  are  given  in  Table  I.  The  RMSEs 
of  the  UKF-E  are  not  shown  in  this  table  for  the  stationary 
platform,  for  the  targets  are  unobservable  in  this  case.  The 

2We  assume  that  the  first  bearing  of  a  target  is  received  from  the  EO/ESM 

sensor,  as  it  is  natural  that  the  non-acoustic  signal  of  a  target  reaches  the 

sensor  platform  earlier  than  the  acoustic  signal. 
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time  (s) 


time  (s) 


Fig.  4.  Maneuvering  platform:  The  estimated  position  RMSE  versus  time 
for  the  target  with  speed  lOm/s. 


Fig.  7.  Stationary  platform:  The  OOSM-AE  estimated  position  RMSE  versus 
time  for  the  targets  with  speeds  of  5m/s  (tgt  1),  lOm/s  (tgt  2),  50m/s  (tgt  3) 
and  lOOm/s  (tgt  4). 


Fig.  5.  Maneuvering  platform:  The  estimated  position  RMSE  versus  time 
for  the  target  with  speed  50m/s. 


Fig.  8.  Stationary  platform:  The  UKF-E  estimated  position  RMSE  versus 
time  for  the  four  targets  with  speeds  of  5m/s  (tgt  1),  lOm/s  (tgt  2),  50m/s  (tgt 
3)  and  lOOm/s  (tgt  4). 


Fig.  6.  Maneuvering  platform:  The  estimated  position  RMSE  versus  time 
for  the  target  with  speed  lOOm/s. 


overall  position  RMSE  for  a  particular  scenario  is  computed 


TABLE  I.  Position  RMSEs  for  all  the  scenarios 


Platform 

Target 

speed 

(m/s) 

Overall  RMSE 

Last  point  RMSE 

OOSM- 

AE 

(m) 

UKF-E 

(m) 

Impro¬ 

vement 

(%) 

OOSM- 

AE 

(m) 

UKF-E 

(m) 

5 

1132.3 

4154.1 

72.7 

114.9 

402.8 

M 

10 

831.6 

4795.7 

85.5 

46.7 

558.9 

50 

806.9 

1530.8 

47.3 

315.5 

755.0 

100 

605.2 

2624.7 

79.9 

171.8 

3975.4 

5 

2724.2 

- 

- 

836.2 

- 

S 

10 

1708.2 

- 

- 

261.3 

- 

50 

755.7 

- 

- 

286.7 

- 

100 

624.2 

- 

- 

445.0 

- 

by 


pos 


RMSE  . 


\ 


,  iV  K 


(55) 


where  i  is  the  run  index,  N  =  100  is  the  number  of  runs, 
K  =  360  is  the  number  of  time  cycles  in  the  scenario,  and 


Pos enW)  =  \/m)-x(^W  +  [y(t^)-y(tl^  (56) 

where  x(t^)  and  are  the  estimated  target  positions,  and 

x(ts^)  and  y(ts^)  are  the  true  target  positions. 
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It  can  be  seen  that  the  OOSM-AE  clearly  outperforms  the 
UKF-E  for  the  maneuvering  platform  scenarios.  The  overall 
accuracy  improvement  in  terms  of  position  RMSE  is  from 
47%  to  86%,  a  significant  improvement.  For  the  slow  moving 
targets  (shown  in  Figs.  3-4),  the  UKF-E  takes  a  longer  time  to 
converge.  The  UKF-E  position  RMSEs  start  to  decrease  at  time 
180s  (after  the  second  turn),  whereas  the  RMSE  reduction  of 
the  OOSM-AE  occurs  around  time  25s,  which  is  much  earlier 
than  for  the  UKF-E.  For  the  fast  moving  targets  (shown  in 
Figs.  5-6),  both  algorithms  converge  fast  at  the  beginning,  but 
the  UKF-E  has  larger  errors  after  a  while. 

For  the  stationary  platform  (Figs.  7-8),  the  OOSM-AE 
provides  reliable  estimation,  whereas  the  UKF-E  diverges  since 
BOT  from  a  single  ESM  sensor  is  not  observable. 

We  also  observe  that  the  OOSM-AE  has  better  performance 
for  the  fast  moving  targets  than  the  slow  moving  targets  in  both 
maneuvering  and  stationary  platform  scenarios.  The  reason 
for  this  is  that  the  slow  moving  targets  have  lower  bearing 
change  rate.  The  information  provided  by  these  slow  changed 
bearings  is  limited  when  they  are  “buried  in  the  noise”,  and 
this  results  in  marginal  observability  and  slow  convergence  at 
the  beginning.  This  effect  is  more  serious  for  the  stationary 
platform  as  its  bearing  change  rate  is  even  smaller  than  for 
the  maneuvering  platform. 

V.  Conclusions 

This  paper  presented  a  new  passive  BOT  approach  through 
fusion  of  an  ESM/EO  and  an  acoustic  sensor  deployed  on 
the  same  sensor  platform.  The  OOSM-AE  algorithm  has 
been  developed  to  estimate  the  target  trajectory  by  utilizing 
the  acoustic  propagation  delay  which  contains  target  range 
information.  This  approach  avoids  the  requirement  for  plat¬ 
form  maneuvers  of  the  conventional  BOT.  Simulation  results 
showed  that  the  OOSM-AE  can  estimate  the  target  trajectory 
effectively  for  the  stationary  platform,  and  provides  significant 
accuracy  improvement  (47%-85%)  over  the  conventional  BOT 
for  the  maneuvering  platform.  This  new  approach  has  signifi¬ 
cant  potential  to  enhance  passive  BOT  capability  significantly. 
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